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Abstract: We re-examine the properties of the Constrained MSSM in hght of updated con- 
straints, paying particular attention to the impact of the recent substantial shift in the Standard 
Model prediction for BR{B Xsj). With the help of a Markov Chain Monte Carlo scanning 
technique, we vary all relevant parameters simultaneously and derive Bayesian posterior probabil- 
ity maps. We find that the case of /j, > remains favored, and that for /j, < it is considerably more 
difficult to find a good global fit to current constraints. In both cases we find a strong preference 
for a focus point region. This leads to improved prospects for detecting neutralino dark matter in 
direct searches, while superpartner searches at the LHC become more problematic, especially when 
/i < 0. In contrast, prospects for exploring the whole mass range of the lightest Higgs boson at 
the Tevatron and the LHC remain very good, which should, along with dark matter searches, allow 
one to gain access to the otherwise experimentally challenging focus point region. An alternative 
measure of the mean quality-of-fit which we also employ implies that present data are not yet con- 
straining enough to draw more definite conclusions. We also comment on the dependence of our 
results on the choice of priors and on some other assumptions. 
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1. Introduction 

Among various possible sets of boundary conditions that one can impose on the multi- 
dimensional parameter space of the effective, low-energy Minimal Supersymmetric Stan- 
dard Model (MSSM), |jl| the by far most popular choice is the so-called Constrained MSSM 
(CMSSM) In the CMSSM, at the GUT scale the soft masses of ah the sleptons, squarks 
and Higgs bosons have a common scalar mass mo, all the gauginos unify at the common 
gaugino mass 1^1/2: ™d so all the tri- linear terms assume a common tri-linear mass pa- 
rameter Aq. In addition, at the electroweak scale one selects tan/?, the ratio of Higgs 
vacuum expectation values and sign(;u), where fi is the Higgs/higgsino mass parameter 
whose square is computed from the conditions of radiative electroweak symmetry breaking 
(EWSB). 

The small number of parameters makes the CMSSM a popular framework for explor- 
ing SUSY phenomenology. Conversely, collider data provides useful constraints on the 
parameter space (PS) of the CMSSM. In the presence of R-parity the lightest neutralino is 
often the lightest supersymmetric particle (LSP). Assuming it to be the dominant compo- 
nent of cold dark matter (DM) in the Universe, allows one to apply the DM relic density 
determination by WMAP and other experiments as a strong constraint on the CMSSM 

PS i. 

Another important constraint on the CMSSM comes from the process b — > 57. An 
approximate agreement of the Standard Model (SM) prediction for BR{B Xs^) with 

^One well-known implementation of the CMSSM is the minimal supergravity model |3|. 
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an experimental determination requires the sum of SUSY contributions, which enter at the 
same 1-loop level, to be strongly suppressed. While the experimental world average has 
for over a year remained at (3.55 it 0.26) x 10~^ j^, the recently re-evaluated SM value of 
BR{B Xgj), as obtained by Misiak et al., in ^, has moved quite substantially from 
(3.60±0.30) X 10"^ down to (3.15 ±0.23) x 10"^. ^ The main shift was caused by including 
new partial NNLO SM contributions, most importantly an approximate evaluation of the 
charm mass effects. The new SM value leads to some discrepancy, at the level of 1.2a, with 
the experimental average. 

From the perspective of SUSY corrections, much more important than this slight 
discrepancy is the fact that the SM central value has now moved from above to below 
the experimental one. In the case of minimal flavor violation, which is applicable to the 
CMSSM, dominant SUSY contributions come from the charged Higgs/top loop, which 
always adds constructively to the SM contribution, and from the chargino/stop loops, 
whose sign is opposite to that of fi. With the previous SM value the 6 ^ 57 branching 
ratio, this was used as an argument for assuming fi to be positive. Indeed, for fi < one 
had to push superpartner masses into the multi-TeV range in order for the chargino/stop 
loop correction to become suppressed, while in the opposite case much smaller masses were 
allowed. However, the recent shift in the SM predictions for BR{B ^s7) makes the 
argument for selecting /i > questionable, and has in fact motivated us to perform this 
analysis. 

Another argument that is often used in favor of /i > is based on a persistent dis- 
crepancy of a^^^* — = (28 it 8.1) x lO^^'^ between the experimental value and the SM 
prediction of the anomalous magnetic moment of the muon {g — 2)^ Taking the nearly 
3.5o" difference as being due to SUSY contributions, ^a^^^"^ = a^^^* — a^, (whose sign is 
the same as that of fi), implies ^ > 0. 

However, such conclusions are based on a somewhat oversimplified treatment of both 
theoretical and experimental uncertainties, which is common practice in fixed-grid scans of 
a SUSY PS. In such studies, a "step- function" approach is usually adopted: regions of the 
PS where contributions to a given observable are within the ztlo" (or some other) range 
around the experimental central value are treated as fully allowed, while those even slightly 
outside are treated as completely ruled out. The same applies to experimental limits, 
e.g., on Higgs or superpartner masses. Instead, it seems more justified to assign varying 
"weights" to different points in a PS, depending on how well, or how poorly, a prediction 
for a given observable matches its experimental determination. Furthermore, in the usual 
approach theoretical errors are typically neglected, and so are residual uncertainties in 
relevant SM parameters, simply for the reason of practicality. 

Recently a more refined procedure has been developed which allows one to overcome 
these shortcomings. It is based on a statistical Bayesian analysis linked with a Markov 
Chain Monte Carlo (MCMC) scanning techniques [ p!0| , 11], and is becoming increasingly 



popular in studying SUSY phenomenology [|12| , |13| , |14| , |15| , |lq , |17|]. The MCMC technique 



further slight decrease to (2.98 ± 0.26) x 10~* after including some additional partial effects due to 
a treatment of a photon energy cut > 1.6 GeV was obtained in ref. Note that the above values do 
depend on the choice of the cut in E-y. 
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allows one to make a thorough scan of a model's full multi-dimensional PS. Additionally 
maps of probability distributions can be drawn both for the model's parameters and also 
for all the observables (and their combinations) included in the analysis. In this approach, 
sharply defined "allowed regions" drawn up in fixed-grid scans are replaced by more infor- 
mative probability distributions. 

The MCMC Bayesian approach to studying properties of "new physics" models, like 
the CMSSM, is superior in the sense of treating the impact of different experimental data 
with their proper weights. It allows to make global scans of the PS and to derive its global 
properties and predictions. When (hopefully) discoveries are eventually made at current 
or future experiments, the approach will provide invaluable in assessing their implications 
for a given theoretical model. 

In this work we apply the MCMC Bayesian formalism to explore the impact on the 
CMSSM's properties from mostly the recent change in the SM value for BR{B Xg'j). 
As we will see, regions of the highest posterior probability, will move rather dramatically 
to the focus point (FP) [l^ region of the CMSSM PS. This in turn will lead to a significant 
shift in prospects for superpartner searches at the LHC (generally for worse) and in direct 
searches for DM neutralinos in the Milky Way (generally for better), while chances of 
finding at the Tevatron will remain good. 

In order to assess the robustness of the results obtained in Bayesian language, following 
our previous work ^ we also apply an alternative measure of a mean quality-of-fit, 
which is similar to a popular x^-measure, which is singles out (possibly limited) ranges of 
parameters that give the best fit to the data. 

We consider both signs of /i. In the probabilistic approach, the case < cannot 
be treated anymore as ruled out, but merely as disfavored, by the {g — 2)^ result. The 
relative weight of this constraint has to be compared with that of other observables in a 
proper statistical way. In a recent similar study of Allanach et al., (although done 
with the old values of BR{B Xg^) and (5a^^^^) fits for both signs of ji were performed. 
It was concluded that the the case < was only marginally disfavored, with the ratio of 
probabilities estimated at P(/x < 0)/P(/i > 0) = 0.07 — 0.16. In our study we also find that 
the case of < gives a worse fit to the data than the opposite sign of /i, although the 
level of preference for > is difficult to quantify. On the other hand, unlike in [^], we 
are not interested in comparing the relative probabilities of the two cases /x < and /U > 0, 
but rather emphasize different implications of each one for various observables of interest. 

In this paper we include, and update when applicable, all relevant experimental con- 
straints from collider direct searches and from rare processes, and also from cosmology on 
the relic abundance of the lightest neutralino O^/i^. We further take into account residual 
error bars in relevant SM parameters. Details of our analysis will be given below. 

We adopt flat priors on the usual CMSSM parameters: mi/2j "^Oj and tan We do 
this primarily for the sake of comparing our results with the literature (in particular with 
the fixed-grid scan approach) where this parametrization is usually assumed. Our specific 
results will accordingly depend on the choice, as we discuss later. 

Implications from the current analysis (assuming /x > 0) for the Higgs bosons have 
already been presented in |16| where we showed that, with our choice of priors, the light- 
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est Higgs boson mass is confined to 115.4 GeV < nifi < 120.4 GeV (95% probability 
interval) and that its couplings to electroweak gauge bosons are very close to those of the 
SM Higgs boson with the same mass. This range should be excluded (at 95% CL) at the 
Tevatron. Here we extend our analysis to to the case /x < 0, reaching similar conclusions 
for at the Tevatron. We also derive most probable ranges of several sparticle masses, of 
the rates of rare bottom quark processes, and of both spin-independent and spin-dependent 
cross sections for dark matter neutralino scattering off nuclei. 

The paper is organized as follows. In Sec. |^ we outline our theoretical setup. In 
Sec. 1^ we present our numerical results for the PS of the CMSSM in terms of the Bayesian 
statistics and of the mean quality-of-fit, and resulting implications for several observables. 
We finish with summary and conclusions in Sec. ^ 

2. The Analysis 

Our procedure based on MCMC scans and Bayesian analysis has been presented in detail 
in |l^]. Here, for completeness, we repeat its main features following an updated presen- 
tation given in [|l6| . 

2.1 Theoretical framework 

In the CMSSM the parameters mi/2, ttlq and Aq, which are specified at the GUT scale 
-^GUT — 2 X 10^^ GeV, serve as boundary conditions for evolving, for a fixed value of 
tan P, the MSSM Renormalization Group Equations (RGEs) down to a low energy scale 
MgusY = ^/^^Wj'W^ (where mj^ denote the masses of the scalar partners of the top quark) , 
chosen so as to minimize higher order loop corrections. At MsusY the (1-loop corrected) 
conditions of electroweak symmetry breaking (EWSB) are imposed and the SUSY spectrum 
is computed at mz- 

We are interested in delineating high probability regions of the CMSSM parameters. 
We consider separately both signs of fi and denote the remaining four free CMSSM param- 
eters by the set 

6* = (mo,mi/2, Ao,tan/3). (2.1) 

As demonstrated in []T2| |l^ , the values of the relevant SM parameters can strongly influence 
some of the CMSSM predictions, and, in contrast to common practice, should not be 
simply kept flxed at their central values. We thus introduce a set ^ of so-called "nuisance 
parameters" of those SM parameters which are relevant to our analysis, 

= {Mt,mb{mt)^, a,m{Mz)^, a,(Mz)^), (2.2) 

where Mt is the pole top quark mass. The other three parameters: m(,{mh)^^ - the bottom 
quark mass evaluated at m^, aem{Mz)^^^ and as{Mz)^^ - respectively the electromag- 
netic and the strong coupling constants evaluated at the Z pole mass Mz - are all computed 
in the MS scheme. 

The set of parameters 6 and ip form an 8-dimensional set m of our "basis parameters" 
m = {9,ip).^ In terms of the basis parameters we compute a number of collider and 

^In |14| we denoted our basis parameters with a symbol rj. 
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cosmological observables, which we caU "derived variables" and which we collectively denote 
by the set ^ = (^^i, ^2, • • •)• The observables, which are listed below, will be used to compare 
CMSSM predictions with a set of experimental data d, which is available either in the form 
of positive measurements or as limits. 

In order to map out high probability regions of the CMSSM, we compute the posterior 
probability density functions (pdf's) p{m\d) for the basis parameters m and for several 
observables. The posterior pdf represents our state of knowledge about the parameters m 
after we have taken the data into consideration. Using Bayes' theorem, the posterior pdf 
is given by 

pHd) = ^«M!!). (2.3) 

p{d) 

On the r.h.s. of eq. (p.3|), the quantity p{d\£,), taken as a function of ^ for fixed data d, is 
called the likelihood (where the dependence of ^{m) is understood). The likelihood supplies 
the information provided by the data and, for the purpose of our analysis, it is constructed 
in Sec. 3.1 of ref. The quantity 7r(m) denotes a prior probability density function 

(hereafter called simply a prior) which encodes our state of knowledge about the values of 
the parameters in m before we see the data. The state of knowledge is then updated to 
the posterior via the likelihood. Finally, the quantity in the denominator is called evidence 
or model likelihood. Here it only serves a normalization constant, independent of m, and 
therefore will be dropped in the following. As in ref. |14| , our posterior pdf's presented 
below will be normalized to their maximum values, and not in such a way as to give a total 
probability of 1. Accordingly we will use the name of a "relative posterior pdf", or simply 
of "relative probability density". 

The Bayesian approach to parameter inference relies on the updating of the prior prob- 
ability to the posterior through the information provided by the data (via the likelihood). 
This requires specification of the prior probabilities for the parameters of the model, that 
in our case are taken to be flat (i.e., constant) over a large range of the CMSSM and SM 
parameters given above. If the data are not strongly constraining, the choice of prior can 
lead to a signiflcant impact through the effect of the "volume" of the parameter space. 
Indeed, as we discussed in ref. imagine the situation that there exist a rather large 
region of the PS where theoretical predictions match the data rather well. In addition, let 
there be a rather small, possibly fined-tuned, region giving very good match of the data. 
The Bayesian posterior probability would give an overwhelming weight to the larger region, 
due to the much larger volume it occupies in parameter space. We notice that this kind of 
situation only arises in the "grey zone" of insufficient data, since of course if the data were 
powerful enough as to rule out such a large region, then the Bayesian posterior probability 
would show this by peaking in correspondence with the best fitting, smaller region. As 



done previously in refs. |14, |16|, we therefore consider also an alternative statistical measure 
of the mean quality-of-fit defined in [p!^], which is much more sensitive to possibly small 
best-fit regions. Below we will compare results obtained using the two measures. 

2.2 Constraints 

We perform a scan over very wide ranges of CMSSM parameters |14, !(:]. In particular 
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SM (nuisance) parameter 


Mean value Uncertainty 
^ a (exper.) 


lei. 


Mt 

as{Mzf^ 
l/aem(Mz)^^ 


171.4 GeV 2.1 GeV 
4.20 GeV 0.07 GeV 
0.1176 0.002 
127.955 0.018 


11 

i 
i 
i 



Table 1: Experimental mean /i and standard deviation a adopted for the likelihood function for 
SM (nuisance) parameters, assumed to be described by a Gaussian distribution. 



we take flat priors on the ranges 50 GeV < mi/2)'"T'0 < 4 TeV (this way including the 
focus point region in this kind of analysis), \Aq\ < 7 TeV and 2 < tan/3 < 62. For the 
SM (nuisance) parameters, we assume flat priors over wide ranges of their values [^4| and 
adopt a Gaussian likelihood with mean and standard deviation as given in table ||. Note 
that, with respect to ref. [14|, we have updated the values of all the constraints.^ 



The experimental values of the collider and cosmological observables that we apply (our 
derived variables) are listed in table with updates where applicable. In our treatment 
of the radiative corrections to the electroweak observables M\y and sin^ ^eff; starting from 



ref. [16 1 we include full two-loop and known higher order SM corrections as computed in 



ref. |^6|], as well as gluonic two- loop MSSM corrections obtained in [27|. We further update 



an experimental constraint from the anomalous magnetic moment of the muon (g — 2)^ 



for which a discrepancy (denoted by 5a^^^^) between measurement and SM predictions 
(based on e^e~ data) persists at the level of 3.5fT We note here, however, that the 
impact of this (still somewhat uncertain) constraint on our findings will be rather limited 
because the corresponding error bar remains relatively large. ^ 

As regards BR{B Xsj), with the central values of SM input parameters as given in 
table 1^, for the new SM prediction we obtain the value of (3.11ib0.21) x 10~^.^ We compute 



SUSY contribution to BR{B Xg^) following the procedure outlined in refs. |2£, 30| which 
were extended in refs. |^, ^ to the case of general flavor mixing. In addition to full leading 
order corrections, we include large tan /3-enhanced terms arising from corrections coming 
from beyond the leading order and further include (subdominant) electroweak corrections. 
Regarding cosmological constraints, we use the determination of the relic abundance of 



^After completing our numerical scans, a new value of the top mass, Mt — 170.9 ± 1.8 GeV, based on 
Tevatron's Run-II 1 fb^^ of data was released |2o|. Including it would not have much impact on our results, 
since the shift in the mean value of Alt is very mild if compared to the standard deviation adopted in this 
paper. 

^Although the different evaluations seem to be converging; e.g., recently (27.6±8.1) x 10^^" was obtained 
in ref. 

^The value of (3.15 ± a2_3) x 10"" originally derived in ref. ||, |l was obtained for slightly different 
values of Mt and a,s{Mz)'^^^ ■ Note that, in treating the error bar we have explicitly taken into account 
the dependence on Mt and as{Mz)^''^ , which in our approach are treated parametrically. This has led 
to a slight reduction of its value. Note also that even though the theoretical error is, strictly speaking, 
not Gaussian, it can still be approximately treated as such as it represents an estimate where a larger 
assumed error of the (dominant) uncertainty due to non-perturbative effects is assigned lower probability - 
M. Misiak, private communication. 
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Observable 



Mw 
siv? 6*2^ 

BR{B Xs-i) X 10^ 



sparticle masses 



Mean value 



Uncertainties 
a (exper.) r (theor.) 



80.392 GeV 

0.23153 

28 

3.55 

17.33 ps-i 
0.104 



29 MeV 
16 X 10"^ 
8.1 
0.26 

0.12 ps-i 
0.009 



15 MeV 
15 X 10"^ 
1 

0.21 

4.8 ps-i 



Limit (95% CL) 



r (theor.) 



< 1.0 X lO"'^ 

> 114.4 GeV (91.0 GeV) 

f{mi) 



14% 
3 GeV 
negligible 



See table 4 in ref. |14| 



ref. 



m 
m 
i 
i 

I 



ref. 



Table 2: Summary of the observables used in the analysis. Upper part: Observables for which 
a positive measurement has been made. (5a^^^^ = p* — denotes the discrepancy between 
the experimental value and the SM prediction of the anomalous magnetic moment of the muon 
{g — 2)^. For central values of the SM input parameters used here, the SM value of BR{B Xgj) 
is 3.11 X 10"'', while the theoretical error of 0.21 x 10"'* includes uncertainties other than the 
parametric dependence on the SM nuisance parameters, especially on Mt and as{Mz)^^ ■ As 
explained in the text, for each quantity we use a likelihood function with mean fi and standard 
deviation s = \/(t^ + t^, where cr is the experimental uncertainty and r represents our estimate 
of the theoretical uncertainty. Lower part: Observables for which only limits currently exist. The 
likelihood function is given in ref. , including in particular a smearing out of experimental errors 
and limits to include an appropriate theoretical uncertainty in the observables. rrih stands for the 
light Higgs mass while = {hZ Z)mssm/ g"^ {hZ Z)sm, where g stands for the Higgs couphng to 
the Z and W gauge boson pairs. 



cold DM based on the 3-year data from WMAP |23| to constrain the relic abundance ri^/i 
of the lightest neutralino which we compute with high precision, including all resonance 



and coannihilation effects, as explained in ref. |14|, and solve the Boltzmann equation 
numerically as in ref. [33|. In order to remain on a conservative side, we impose the 
following Gaussian distribution 



n,^h^ = 0.104 zb ^(0.009)2 + (0.1 Qy.h'^)'^ = 0.104 zb 0.009 yi + 1-335 (17^/i2/0.104)2(2.4) 

Note that our estimated theoretical uncertainty is of the same order as the uncertainty 
from current cosmological determinations of f^cDM^^- 

We further include in our likelihood function an improved 95% CL limit on BR{Bs — > 
H^H^) and a recent value of Bg-Bs mixing, AM^^, which has recently been precisely 
measured at the Tevatron by the CDF Collaboration p^] . In both cases we use expressions 



from ref. |32] which include dominant large tan /^-enhanced beyond-LO SUSY contributions 
from Higgs penguin diagrams. Unfortunately, theoretical uncertainties, especially in lattice 
evaluations of Jb^ are still very large (as reflected in table |2| in the estimated theoretical 
error for AM^^), which makes the impact of this precise measurement on constraining the 
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CMSSM parameter space rather limited.'' 

For the quantities for which positive measurements have been made (as hsted in the 
upper part of table we assume a Gaussian likelihood function with a variance given by 
the sum of the theoretical and experimental variances, as motivated by eq. (3.3) in ref. 
For the observables for which only lower or upper limits are available (as listed in the bottom 
part of table 0) we use a smoothed-out version of the likelihood function that accounts for 



the theoretical error in the computation of the observable, see eq. (3.5) and fig. 1 in |14|. 
In applying the Higgs boson lower mass bounds from LEP-II we take into account 
its dependence on its coupling to the Z boson pairs = g'^ {hZ Z)mssm/ d'^ {hZ Z)sm, as 
described in detail in ref. |16|. When (^^ ~ 1, the LEP-II lower bound of 114.4 GeV 



(95% CL) is applicable. For arbitrary values of C/n we apply the LEP-II 95% CL 
bounds on m/j and m^, which we translate into the corresponding 95% CL bound in the 
(mhXh) plaiie- We then add a theoretical uncertainty T{mh) = 3 GeV, following eq. (3.5) 
in ref. ||l^. However, a posteriori we find ~ 1 which means that the CMSSM light Higgs 
boson is invariably SM-like. This procedure results in a conservative likelihood function 
for nifi, which does not simply cut away points below the 95% CL limit of LEP-II, but 
instead assigns to them a lower probability that gradually goes to zero for lower masses. 

Finally, points that do not fulfil the conditions of radiative EWSB and/or give non- 
physical (tachyonic) solutions are discarded. We adopt the same convergence and mixing 



criteria as described in appendix A2 of ref. |14], while our sampling procedure is described 



in appendix Al of ref. []14| . We have the total of = 10 MC chains, with a merged 
number of samples 2 x 10^, and an acceptance rate of about 1.5%. We adopt the Gelman 
&: Rubin mixing criterion, with the inter-chain variance divided by the intra-chain variance 
(the R — 1 parameter) being less then 0.1 along all directions in parameter space. More 
details of our numerical MCMC scan can be found in [14|. 



3. Results 

We will now explore the implications of the above constraints on the CMSSM parameters, 
paying some more attention to the impact of BR{B Xg^). We will compare the posterior 
probability distributions of the Bayesian language with the ranges favored by the mean 
quality-of-fit. Next, we will discuss implications for Higgs and superpartner masses and for 
direct detection of the neutralino dark matter. In computing the Higgs (and SUSY) mass 



spectrum we employ the code SOFTSUSY v2.08 |l35 |. 



3.1 Implications for the CMSSM parameters 

We first show in fig. |l| the 2-dim relative probability density functions in the planes spanned 
by the CMSSM parameters: mi/2, tuq, tan/3, Aq, and assuming /u > 0, while in fig. ^ the 
same is shown for fi < 0. In each panel all other basis parameters have been marginalized 
over. Redder (darker) regions correspond to higher probability density. Inner and outer 



^On the other hand, in the MSSM with general flavor mixing, even with the current theoretical uncer- 
tainties, the bound from AMb^ is in many cases much more constraining than from other rare processes 
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Figure 1: The 2-dim relative probability density functions in the planes spanned by the CMSSM 
parameters: mi/2, fno, Aq and tan/? for /i > 0. The pdf's are normalized to unity at their peak. 
The inner (outer) blue solid contours delimit regions encompassing 68% and 95% of the total 
probability, respectively. All other basis parameters, both CMSSM and SM ones, in each plane 
have been marginalized over. This figure should be compared with figure 2 in ref. 



blue (dark) solid contours delimit regions of 68% and 95% of the total probability, respec- 
tively, and remain well within the assumed priors, except for mg. In all the 2-dim plots, 
the MC samples have been divided into 70 x 70 bins, with a mild smoothing across adja- 
cent bins to improve the quality of the presentation (this has not impact on our statistical 
conclusions). Jagged contours are a result of a finite resolution of the MC chains. 

In the case of ;U > (fig. |l|) we can see a strong preference for large mo > 1 TeV. On 
the other hand, the peak of probability for mi/2 is around 0.5 TeV, although the 68% range 
of total probability is rather wide, increases with mo and exceeds 1.5 TeV for niQ ~ 4 TeV. 
Additionally, at smaller mo < 1 TeV there are a few confined 68% total probability regions. 

The strong preference for large mo ^ ^1/2 is primarily the result of the sizable shift 
in the SM value of BR{B — > Xgj), as can be seen by comparing fig. || with fig. 2 in ref. |jl^ 
(or fig. 8 of ref. [^]) where the previous value of BR{B ^s7) has been used. (While 
the other CMSSM parameters also experience some shift in their most probable values, it 
is not as dramatic as that of mo towards larger values.) The underlying reason is that, 
at fairly small mi/2 the charged Higgs mass remains relatively Hght, in the few hundred 
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GeV range, and, via a loop exchange with the top quark, it adds substantiahy to the SM 
value of BR[B ^sT), towards the experimentally allowed range. (In fact, for m^^- ~ 
650 GeV, the contribution is sufficient to fill the gap between the SM and the experimental 
central values of BR{B — > Xsj) Q.) At smaller mi/2 and/or mg, the (negative, for 
/i > 0) chargino-stop contribution is too large and needs to be compensated by the H~-top 
contribution. In fact, we do find some small "islands" of 68% total probability at mi^2 ^ 
1 TeV and mo < 1.3 TeV (in particular, notably, an interesting case of — 0-5 TeV 
and rriQ ~ 0.2 TeV) but the bulk of high probability region corresponds to tuq > 1 TeV. 

At 95% total probability level the available parameter space opens widens considerably, 
and also some new features arise. In particular, at mi/2 — 0-2 TeV we can see a narrow 
high-probability funnel induced by the light Higgs boson resonance [12, 14]. Also, tan/3 
becomes less confined to its most preferred range of large values between some 50 and 60, 
while ^0 remains on a positive side. 

In the case of fx < (fig. ^), one can see a strong preference for even larger mo. Also, 
the 68% total probability region of mi/2 shifts towards larger values, although still remains 
basically below 2 TeV. (Although we again find an interesting isolated high probability 
region at m^/2 — ^■'^'^ TeV and mo — 0.2 TeV.) This shift towards larger mo and/or 
mi/2 is again caused mostly by the BR{B — > ^^7) constraint. At large mg (and not too 



-10- 




Roszkowski, Ruiz & Trotta (2Q07) 




-4 -2 2 4 6 



Ao (TeV) 



Roszkowski, Ruiz & Trotta (2007) 




20 40 60 

tan|8 



Figure 3: The 1-dim relative probability densities for the CMSSM parameters mo, nT.i/2, and 
tan /3. All other parameters have been marginalized over. Solid blue (dashed red) curves correspond 
to /i > (/i < 0) . This figure should be compared with figure 4 (black solid lines) in ref . . 

large mi/2) the charged Higgs mass decreases and its contribution tends to be on a high 
side, while the chargino-stop one becomes too small. A similar effect is observed at large 
mi/2 ~ 1-5 TeV and mo > 1 TeV where the chargino and stop masses become too large to 
contribute much as well. 

In fig. ^ the 1-dim marginalized probability distributions for the CMSSM parameters 
are compared for both > and /x < 0. In each panel all the other CMSSM parameters and 
all SM (nuisance) parameters have been marginalized over. It is clear that non-negligible 
probability ranges of the CMSSM parameters, other than mo, are confined well within their 
assumed priors. Again, we can see strong preference for large mo ^ 1 TeV (even stronger 
for /X < than for > 0). Larger values of ^-re also favored for /x < although, 
for both signs, this parameter is well confined within 2 TeV. The preferred range of Aq is 
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fairly uncorrelated with the other parameters |14], and it is symmetrically peaked (basically 
independently of the sign of /i) at some 1.5 TeV. This value is however basically twice as 
large as for the previous SM value of BR{B — > Xg^) [^]. On the other hand, tan/3 is 
well-peaked at some 53 for /i > and some 48 for < 0. In both cases, there remains a 
sizable tail of much smaller values which remain allowed at large mo- 

What most strongly contributes to confining '^^11 within its prior (for both signs 
of /i) is the relic abundance fi^/i-^ which becomes too large. On the one hand, at large 
tan/3 it becomes easier to satisfy the constraint from ^-xh? due to the increased role of the 
neutralino annihilation via the pseudoscalar Higgs effect and/or the coannihilation effect. 
On the other, as explained in ref. [0, as tan/3 becomes very large, > 60 for /u > 
(> 50 for /i < 0), it becomes very difficult to find self-consistent solutions of the RGE's. 

The feature that very large values of mo (the FP region), up to the assumed prior of 
4 TeV, remain allowed (actually, even preferred), is unfortunate but is a consequence of 
the fact that current data are not constraining enough. Normally, as all the superpartner 
masses (including the LSP) increase, bino-like neutralino annihilation in the early Universe 
becomes suppressed and it becomes harder to satisfy the WMAP constraint on Q.^h? . This 
is why mi 12 is well confined within some 2 TeV, as described above. Unfortunately, in 
the FP region the behavior of Q^h^ is much more sensitive to input parameters. We 
illustrate this in fig. ^ where we plot Q^h? vs. mo for /_i > and a choice of the other 
CMSSM parameters close to their highest probability values (fig. |l]). Clearly, as Mt is 
varied within la around its central value (cv), Vl-^h? changes quite substantially. In the 
presented example, by fixing Mt at its central value [^6| (as it is normally done in fixed- 
grid scans) one would find no cosmologically allowed mo. On the other hand, by reducing 
(increasing) Mt by la we can find one (two widely disconnected ) narrow region(s) of mo 
where fi^/i^ in the WMAP range. Furthermore, in a probabilistic approach, even the case 
at the central value of Mt is not excluded but only less favored by Q.^h? . 

Another feature that is evident in fig. ^is that, as Mt is varied by la around its central 
value, the range of mo where self-consistent solutions of the RGE's and the conditions of 
EWSB can be found changes by as much as a factor of two. It is therefore clear that, if 
one includes the FP region, it is basically impossible to locate the cosmologically favored 
range of mo. In particular, it would be misleading to simply fix the top mass at its central 
value (and likewise with the bottom mass at large tan/3). 

The rather special properties of the FP region are to a large extent related to the 
behavior of the parameter /i. As mo increases (along, say, fixed mi/2); /^^ (which is deter- 
mined the the conditions of EWSB) decreases rather quickly, thus increasing the higgsino 
component of the neutralino. This in turn reduces Q^^h? to an acceptable range for some 
narrow range of mo. At slightly larger mo, /i^ drops below zero, thus delimiting the zone 
where consistent, non-tachyonic solutions can be found. Thus generally in the FP region 

is rather small relative to mo- We can see this feature in fig. |5| where we plot the 1- 
dim relative probability density for the parameter. We can see a clear peak in the few 
hundred GeV region for both signs of ^. In this region ~ m^/2/3- Additionally, for 

< there is a second well-pronounced peak around some 1.4 TeV which corresponds to 
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Figure 4: An illustration of the sensitivity of the cosmologically favored range of toq to the top 
mass Mt in the focus point region. Other CMSSM parameters have been fixed close to their most 
preferred values for /.t > (see fig. |l|). We show ^x^"^ (blue solid line) for the central value (cv) of Mt 
(middle panel) and for the values decreased and increased by la (left and right panel). The parallel 
dotted lines denote the very narrow 2a range, as determined from the 3-year data of WMAP p3| . 
The long-dashed line denotes the combined theoretical plus experimental error, as described in the 
text. Note that, as Mt is varied between its cv minus \a to its cv plus la, the range of mo, for 
which consistent solutions can be found, roughly doubles. 
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Figure 5: The 1-dim relative probability density of All other parameters have been marginal- 
ized over. Dashed blue (dotted red) curves correspond to > (/l( < 0). 

the band of higher relative probability around 771-1/2 — 1-5 TeV and mo > 1 TeV in fig. |2|, 
outside of the FP region. 

In order to summarize the above discussion, in table ^ we give the 68% and 95% total 
probability ranges of the CMSSM parameters for both signs of /i. 

We emphasize that the above results do not at all imply that it is equally easy to fit 
all the data for both signs of ^. (In figs. ||- ^ the posterior probabilities are normalized 
relative to the respective highest values.) Actually, for negative jx the fit is considerably 



-13- 



Parameter 


/i < U 

68% region 95% region 


/i > 

68% region 95% region 


niQ (TeV) 
mi/2 (TeV) 
Ao (TeV) 
tan P 


< 3.51 < 3.93 
(0.57,1.55) (0.87,1.86) 
(-0.19,3.42) (-1.79,5.49) 
(18.6,46.08) (7.51,48.9) 


< 3.1 < 3.87 
(0.4,1.16) (0.145,1.6) 
(0.11,2.94) (-1.86,4.84) 
(26.38,54.18) (11.17,56.78) 



Table 3: CMSSM parameter ranges corresponding to 68% and 95% of posterior probability (with 
all other parameters marginalized over) for both signs of /i. 



poorer. In fact, we find that the best-fit foi' the ^ > case is 6.3, while for the ^ < 
case it more than doubles to 14.4. It is difficult to attach a precise statistical significance 
to this result, as clearly the distributional properties of the parameter space are far from 
being Gaussian (and hence the is not chi-square distributed). A proper evaluation of the 
significance of this difference in the goodness-of-fit would require Monte Carlo simulations 
of the measurements for both the > and /x < cases, which is beyond the scope of this 
work. However, this is an indication of the fact that the /x < case is at greater tension 
with the data than /i > 0, although it cannot be conclusively ruled out yet. 

A fully Bayesian approach would consider computing the Bayes factor among the two 
possibilities for /i, along the lines of what has been done in ref. [|l5|. This procedure 
is, however, computationally demanding, and the result is potentially strongly dependent 



on volume effects deriving from the choice of priors (see, eg., ref. |37]). An interesting 
alternative is to use the procedure outlined in ref. [^] which employs Bayesian calibrated p- 
values to obtain an upper limit on the Bayes factor regardless of the prior for the alternative 
hypothesis. In the present case, the application of this procedure would require Monte 
Carlo simulation to obtain the p-value corresponding to the observed x-SQuare difference. 
However, if we take again the result of Ax^ = 7 at face value, assuming that it is indeed 
X-square distributed (which is probably a very poor approximation, as argued above), then 
the corresponding upper limit on the Bayes factor is 10 : 1. (This rough estimate is actually 
in surprising good agreement with the result found in ref. [^] using numerical integrations) 
This would mean that the minimum probability of ^ < is 0.1, which certainly does not 
constitute strong evidence against /i < 0. The above considerations only highlight the 
difficulty of translating our result into a precise statement about the relative probability 
of /i > vs /X < 0. 

A related (although somewhat different) issue is to identify regions of the CMSSM PS 
where the fit to the data is much better than elsewhere. As stated above, if such regions 
occupy a small volume of parameter space (given our choice of priors), the posterior pdf 
will consequently weight them down, although they might exhibit a higher goodness-of-fit. 
To address this point, we consider an alternative measure of the mean quality-of-fit, which 
is defined as the average of the effective under the posterior distribution, i.e.. 



{X"") = J dm[-2\np{d\m)]p{m\d), (3.1) 
which is a quantity that is largely insensitive to the choice of priors (as long as the best- 
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fitting points are explored by the MCMC scan). Its distribution for tlie CMSSM parameters 
is plotted in figs. ^ and ^ for // > and /i < 0, respectively, in each case they are normalized 
to the respective best-fit value. We can see that, for /i > 0, there indeed exists at least one 
well-localized region around 771-1/2 — 0-4 TeV and mo — 1-5 TeV (not far from the location 
of the highest relative pdf, and in any case within the 68% total probability contour), with 
another one at somewhat smaller values of both mi/2 and niQ. (Below we will show however 
that such best-fit regions may be in conflict with dark matter search limits, which have 
not been applied as constraints at this stage.) On the other hand, at < it is generally 
more difficult to find a good fit to the data, as indicated by the larger value of the given 
above. We also notice that many of the best-fitting regions for the case < in fig. ^ lie 
outside the 95% posterior probability contour, further indicating a strong tension with the 
data. 
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Figure 6: The mean quality-of-fit in the planes spanned by the CMSSM parameters: mi/2, Tip, 
Aq and tan/3 for /x > 0. For comparison, (blue solid) total posterior probability contours of 68% 
and 95% from fig. [| have been added. This figure should be compared with figure 11 in ref. jl^ . 

Different experimental observables may constrain or favor different regions of the 
CMSSM parameters; they may "pull" in different directions. We display this in fig. ^ 
where we plot the 1-dim posterior pdf's for several variables for both signs of \x. For com- 
parison, we also plot the corresponding Gaussian likelihood functions representing the data 
used in the fit. If there was no tension among different observables then, in the absence 
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Figure 7: The same as in fig. g but for /i < 0. For comparison, (blue solid) total posterior 
probability contours of 68% and 95% from fig. || have been added. 



of strong correlations among them, the relative probability curves should overlap with the 
data. This is basically the case for Mt, as{Mz)^^^ and Q.yjT? . (In the last case the slightly 
skewed shape of the pdf 's is a result of our treatment of the theoretical uncertainty which 
is larger for larger Q.-^h?.) On the other hand, the electroweak observables my/ and sin^ 6'eff 
show some pull away from their expected values, in general agreement with ref. [^] where 
the two variables were computed with a similar precision. On the other hand, the tension 
is insufficient to provide convincing preference for low MgusY, in apparent contrast to the 
findings of ref. [^] . 

The biggest tension between best-fit values and experiment is displayed, unsurprisingly, 
in the SUSY contribution to the anomalous magnetic moment of the muon. For both signs 
of n the peaks of the relative probability are far below the central experimental value (about 
3.2(7 for ;U > versus about 3.7cr for /i < 0), and close to each other We conclude 

that it is not justified to use this sole observable to select the positive sign of /i - one has 
to perform a global fit in all of the variables and judge the two cases by this criterion. 
Furthermore, the new BR{B Xsj) actually seems to agree with the data slightly better 
for /i < than for the other sign. Generally, for fi > the total BR{B — > Xs^) remains 
peaked around the SM central value, while for < it is somewhat above it. Finally, 
BR{Bs — > and AMb^ are peaked at their SM values, somewhat more so than 
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Figure 8: The l-dim relative probability density functions for several variables: Mt, as[MzY'^^ j 
n^h"^, Amw = mw - 80.392 GeV, Asin^e'eff = sin^ 6'eff - 0.23153, BR(B Xsi), Saf'^^ , 
BR{Bs iX^iJL~), AMb,- In each panel the dashed blue (dotted red) curves correspond to > 
(/i < 0) while the black solid line represents the data as encoded in the likelihood function. This 
figure should be compared with figure 10 in rcf. 11. 



before ||T4|. 

3.2 Implications for collider searches 

In our Bayesian formalisms, high probability ranges of the CMSSM parameters can easily 
be translated into analogous ranges for Higgs and superpartner masses, and for other 
observables, including indirect processes and dark matter detection cross sections. We 
discuss them in turn below. 

We start with the Higgs bosons. In fig. |^ we plot the Bayesian relative probability 
density distributions of the mass of the lightest Higgs boson and of the charged Higgs 
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Figure 9: The 1-dim relative probability density functions for the masses of (left panel) and 
(right panel). Dashed blue (dotted red) curves correspond to > (/i < 0). Analogous and 
more detailed figures for the case /i > were presented in fig. 7 of ref. Ilfl]. 



boson for both signs of fi.^ (The other Higgs bosons are basically degenerate in mass 
with H^.) In both cases we can see a clear peak in rrih close to 120 GeV, and a sharp 
drop-off for larger values of the mass. The other Higgs are typically considerably heavier 
for negative fi than for the other choice. This may provide one way of an experimental 
determination of the sign of ^. The 68% and 95% total probability ranges of the Higgs 
masses are given in table |^ below. We should also note that the alternative measure of the 
mean quality-of-fit favors lower ranges of the masses of all the Higgs bosons in the case 
of /.i > For the case fi < 0, the mean quality-of-fit distribution of rrih is roughly 

similar to that of the pdf but with the peak shifted to the right by about 1 GeV (plus 
some moderate preference for smaller values, around the current LEP-H limit), while the 
masses of all the other Higgs bosons are preferably very heavy, in the TeV regime. This 
is a reflection of the absence of regions giving good fit to the data in the case of /i < 
(compare fig. |^. 

In ref. [|l6| we have investigated in detail light Higgs masses and couplings for the case 
of /i > 0. In particular we showed that, throughout the whole CMSSM parameter space, 
the couplings of the lightest Higgs boson to the gauge bosons Z and W are very close 
to those of the SM Higgs boson with the same mass, while its couplings to bottom quark 
and tau lepton pairs show some variation. We concluded that, at the Tevatron, with about 
2fb~^ of integrated luminosity per experiment (already on tape), it should be possible to 
set a 95% CL exclusion limit for the whole 95% posterior probability range of nih- Based 
on fig. ^ and table ^ we extend this conclusion to the case of /i < 0. On the optimistic side, 
should a Higgs signal be found, in order to be able to claim a 3a evidence, at least about 
4fb~^ will be needed, independently of the sign of fi. The Tevatron's ultimate goal is to 

*The case of /i > was already presented in fig. 7 of ref. mM. We include it here for comparison. 
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collect about 8fb per experiment. 

One should remember that the above conclusions do depend on the assumed prior 
range of niQ < 4 TeV, as well as on the choice of adopting flat priors in the CMSSM 



variables of eq. 2.1. For instance, adopting a much more generous upper limit mo < 8 TeV 
would lead to changing the ranges for > to roughly 120.4 GeV < mh < 124.4 GeV 
(68% CL) and 115.4 GeV < m,, < 125.6 GeV (95% CL), the latter of which could be 
excluded at 95% CL with about 3fb~^ of integrated luminosity per experiment [^]. Still, 
should no Higgs signal be found at the Tevatron, large ranges of tuq will become excluded 
at high CL, with the specific value depending on the accumulated luminosity. 
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Figure 10: As in fig. but for the masses of several representative superpartners. This figure 
should be compared with figure 5 (black solid lines) in ref. [l^ . 

Turning next to superpartners, we show in fig. |l^ the relative pdf's of the masses of 
several of them, while in table ^ we give the corresponding 68% CL and 95% CL ranges. 
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Table 4: Higgs boson and selected superpartner mass ranges (in TeV) containing 68% and 95% 
of posterior probability (with all other parameters marginalized over) for both signs of /i. Masses 
above 1 TeV have been rounded up to 1 significant digit. 



The blue dashed curves are for /x > and the red dotted ones for < 0. Firstly, for 
/X < all the scalar superpartners are considerably heavier than for the other sign of /x, as 
expected based on the discussion of most probable ranges of mi/2 and especially tjiq. In 
fact, ii fi < then all the sleptons and squarks (whose masses, except for the 3rd generation, 
are at least as large as mg) may be beyond the reach of the LHC. For > there is a 
good chance of seeing the gluino (assuming the LHC reach of some 2.7 — 3 TeV) and a 
reasonable chance of seeing some squarks and sleptons. Unfortunately, these prospects are 
considerably less optimistic than what we found in fig. 5 and table 6 of ref. |14|, where the 



previous SM value of BR{B Xgj) was used. A dedicated analysis would be required to 
derive more detailed conclusions. 

3.3 Implications for direct detection of Dark Matter 

We will now examine implications for direct detection of the lightest neutralino assumed 
to be the DM in the Universe, via its elastic scatterings with targets in underground detec- 
tors Q. We will consider both spin-independent (SI) and spin-dependent (SD) interactions. 
The underlying formalism for both types of interactions can be found in several sources. 
(See, e.g., |3^, |4^, 0.) In this analysis we use the expressions and inputs as presented 
in ref. |41]. We only note here that the SI interactions cross section 0"^^ of a WIMP scat- 



tering off a proton in a target nucleus is the same as that of a neutron and that the total SI 
interactions cross section of the nucleus is proportional to ^ times the square of the mass 
number. In contrast, for the SD interactions, the cross section for a WIMP scattering off a 
proton, cTp^, does not necessarily have to be the same as the one from a neutron |E1, 43]. 



In fig. 11 we show the Bayesian posterior relative probability distribution in the usual 



plane of 0"^^ and the DM neutralino mass for /x < (left panel) and /x > (right 
panel). Starting with /x > 0, we can see a big concentration of probability density at rather 
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Figure 11: The 2-dim relative probability density for dp vs. the neutralino mass for /i < 
(left panel) and /i > (right panel). The inner (outer) solid contours delimit the regions of 68% 
and 95% total probability, respectively. Some current experimental upper bounds are also shown. 
The right panel should be compared with figure f3 (top) in ref. Q. 
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Figure 12: The mean quality-of-fit for a^^ vs. the neutralino mass for /i < (left panel) and 
/J, > (right panel). The inner (outer) solid contours delimit the regions of posterior 68% and 95% 
total probability, respectively (compare fig. |l^). Some current experimental upper bounds are also 
shown, but have not been included as constraints in the likelihood function. The right panel should 
be compared with figure 13 (bottom) in ref. 



high values of cXp^ ~ 10~^ pb, characteristic of the FP region of large mo [44|, which is 
favored by the current theoretical evaluation of BR{B ^sj), as we have seen above. In 
the FP region, the neutralino, while remaining predominantly bino-Iike, receives a sizable 
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higgsino component, which strengthens the dominant Higgs-exchange contribution to ■ 
In addition, there is a more weh- known branch, with cjp^ decreasing with m^, which comes 
from the (now somewhat disfavored) region of rrii /2 , mo ^ 1 TeV where the rehc abundance 
O^/i^ of the neutrahnos is reduced to agree with WMAP and other determinations by a 
pseudoscalar Higgs resonance in their pair annihilation and/or by their coannihilations 
with sleptons. In order to appreciate the change in the CMSSM predictions for fjp^, the 
right panel should be compared with the top panel of fig. 13 in ref. |14] where a previous 
value of the SM prediction for BR{B — > Xgj) was used. 

The left panel of fig. ^ (the case of ;U < 0) also shows a high preference for ~ 
10^*^ pb, which corresponds to the FP region at multi-TeV mo, as for the other sign of fx. In 
addition, we find another rather large 68% total probability region at extremely low values 
of below 10~^'' pb, which corresponds to the higher probability region of mi/2 ~ 1.5 TeV 
in the (mi/2;"^o) plane. Such tiny ranges of are a result of cancellation between the 
Higgs-exchange contribution to up- and down-type quarks 4C]. 

In both panels of fig. 11 we have marked some of the current direct experimental upper 
limits [45, assuming a default value of 0.3 GeV/cm^ for the local DM density. It 

is encouraging that experiments, notably XENON-10 with its very new limit, are already 
probing some portions of the CMSSM PS for /i > 0. CDMS-II is currently taking data 
and is expecting to improve its limit to a similar level of sensitivity. Clearly, a further 
improvement by about an order of magnitude will constitute a critical leap as it will 
hopefully allow one to reach down to the heart of the SI interactions cross sections favored 
in the CMSSM. Future one-tonne detectors are expected to reach down to cjp^ > 10"^'^ pb, 
thus probing most of the favored parameter space of the CMSSM, at least for the more 
favored case of /i > 0. We note that the probability of 0"^^ > 10^^*^ pb in fig. 11 is 98.4% 







for /i > and 62.5% for ^ < 0. 

It is worth re-emphasizing that, in addition to light Higgs searches at the Tevatron 
and the LHC, direct detection DM searches will provide an alternative experimental probe 
of the FP region at multi-TeV mo, in which the squarks will be too heavy to be produced 
at the LHC, although the gluino mass should remain mostly (partly) accessible for > 
ifi < 0). 



For comparison, fig. |12| shows the ranges favored by the alternative measure of the mean 
quality-of-fit. Starting with the case of /.i > (the right panel), we can see a handful of 
small regions of a rather large 0"^^, above a few times 10""^ pb, which are already in conflict 
with the current limit from the XENON-10. These cases corresponds to the "islands" of 
good fit to the data that we have already seen in fig. ^ where they are all concentrated 
in the region of mi^2 ^ 0-5 TeV and mo < 1.5 TeV. With a modest improvement of 
sensitivity, DM search experiments will be able to probe the entire region favored by the 
mean quality-of-fit in the case of ^ > 0. (Notice that neither the posterior pdf nor the 
mean quality-of-fit give much preference for very small .) In contrast, for fi < (the 
left panel of fig. ^) there is hardly any region giving good quality fits. This is a reflection 
of what we have already seen in the (mi/2;^o) plane in fig. |^. 

Turning next to SD interactions, in fig. 13 we present the relative probability density for 
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Figure 13: The 2-dim relative probability density of ap vs. the neutralino mass for /i < 
(left panel) and /i > (right panel). The inner (outer) solid contours delimit the regions of 68% 
and 95% total probability, respectively. Some current experimental 90% CL upper bounds are also 
shown. Analogous plots for are basically identical but the experimental limits from direct 
detection are weaker by nearly two orders of magnitude. 



neutralino-proton scattering cross section cjp^ versus for /i < (left panel) and > 
(right panel) . In the FP region the increased higgsino component of the neutralino in this 
case leads to a larger coupling to the dominant Z-boson exchange. This is reflected in the 
figure where the highest probability regions are, for both signs of /U, concentrated around 
(Tp ^ ~ 10"'* pb. For fi > there is an additional higher probability region which is visible 
in the right panel, and which corresponds to the Higgs resonance and coannihilation region 
mentioned above. In the case of ^ < instead, an additional higher probability region is 
visible in the left panel at a'^^ < 10~^pb and ~ 0.7 TeV. It corresponds to the region 
of 1711/2 ~ 1-5 TeV in the {mi/2,mo) plane (fig. |2|). 

The current experimental upper limits |4^, |5^, 51 1 from direct searches, assuming 
a default value of 0.3 GeV/cm^ for the local DM density, as well as an indirect limit from 
neutralino annihilations into neutrinos in the Sun, the Earth or the Galactic center |B2|, 



which have also been shown in fig. |13|, still remain a few order of magnitude above the 
predictions of the CMSSM. On the other hand, experimental sensitivity has undergone 
steady progress also in the case of SD interactions. Eventually, it will be important to 
reach down below the level of < 10~^ pb, which would allow one for an independent 
cross-check of CMSSM predictions for dark matter. 

In table |5| we have listed the ranges of both cr^^ and a'p^ containing 68% and 95% 
of posterior probability (with all other parameters marginalized over) for both signs of fi. 
The ranges cover the whole allowed range of and provide supplementary information 
to what one can read out of figs. ^ and 13. 

As mentioned above, in the SD interaction case the cross section of WIMP scattering 
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SDiri-iTirlpDPTirlpnt, cv 
68% 


nss sprtion rr^^ ( viW\ 
95% 




(2.9 X 10~^2 81 X 10"^) 


(1.6 X 10-1^, 1.4 X 10-8) 




(2.8 X 10-1°, 3.9 X 10-8) 


(7.2 X 10-i\2.5 X lO-'^) 




Spin-dependent cro 
68% 


3S section a^^ (pb) 
95% 




(5.3 X 10-^9.6 X 10-5) 


(1.7 X 10-8,4.4 X 10-^) 


fi>0 


(2.0 X 10-^8.9 X 10-5) 


(4.2 X 10-8,4.9 X 10-^) 



Table 5: Direct detection of dark matter: ranges of spin-independent and spin-dependent cross sec- 
tion corresponding to 68% and 95% of posterior probability (with all other parameters marginalized 
over) for both signs of /j,. 



from a neutron can in general be very different from the one from a proton. However, in 
our scan of the CMSSM, we have found its relatively probability distribution versus 
to be very close to the case with the proton, and thus do not show it here. As regards 
the mean quality-of-fit (not shown for cr^^), at // > the best-fit regions in fig. ^ give a 
rather small value of ~ IQ-^ pb, while at /u < no range of a^^ is favored by the mean 
quality-of-fit criterion. 
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Figure 14: The 2-dim relative probability density of vs. u^/ for /i < (left panel) and /i > 
(right panel). The inner (outer) solid contours delimit the regions of 68% and 95% total probability, 
respectively. All other parameters have been marginalized over. 



3.4 Correlations among observables 

The probabilistic approach employed in this analysis makes it very easy to examine various 
possible global correlations among different observables in the CMSSM. For example, in 
fig. 14 we present the 2-dim pdf of a^^ versus a^^ for ;U < (left panel) and /i > 
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Figure 15: The 2-dim relative probability density for pairs of selected variables for /i > 0. The 
inner (outer) solid contours delimit the regions of 68% and 95% total probability, respectively. All 
other parameters have been marginalized over. This figure should be compared with figure 14 in 
ref. [pl. 



(right panel). For both signs of ^ we can see a big concentration of probability density on 
cTp^ ~ 10^"® pb and a^^ ~ 10~^pb, which is a reflection of their behavior in fig. 11) and 



fig. |13|, respectively - an effect of the FP region. For > we can see an additional feature 
of a positive correlation which is due to the contribution from the Higgs resonance and/or 
coannihilation effects. 



More correlations are displayed in fig. 15 for the case of /i > 0. Again, the effect 



of the FP region is overwhelming. In all the observables, other than the dark matter SI 
cross section, SUSY effects likely to be tiny ~ the probability density is clearly peaked at 
the respective SM values. For fi < (not displayed) the concentration is typically even 
stronger. Fig. 15 should be compared with fig. 14 of ref. |14| in order to appreciate the 



change. Some correlations which were quite well pronounced in that figure (eg. in a^^ 
versus BR{Bs fi^ fi~)) are now barely visible. 

4. Summary and conclusions 

We have applied the highly efficient MCMC scanning method to explore the parameter 
space of the CMSSM and outlined the regions favored by the Bayesian posterior probability 
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and, for comparison, by the mean mean quality-of-fit. We assumed flat priors in the 
usual CMSSM parameters, applied and updated all relevant experimental constraints from 
colliders and cosmological dark matter abundance, while paying particular attention to the 
impact of the recent change in the SM value of BR{B Xg^f). We examined both signs 
of /i. For both choices, we found strong preference for the focus point region with large 
rriQ in the few TeV range (for ^ < actually saturating at the assumed prior boundary of 
4 TeV), and not as large mij2 < 2 TeV. For comparison, the mean quality-of-fit measure 
selected a small number of isolated regions giving good fit to the data for > but not 
for /i < 0. It appears that the choice /i < is more at odds with the data, as refiected by 
its worse quality-of-fit. 

We then examined ensuing implications for Higgs and superpartner searches and for 
direct detection of dark matter. Prospects for the Tevatron of excluding the whole 95% CL 
range of m/j remain very good, or else there is very good hope to see at least some evidence 
of a signal with the expected final integrated luminosity. Scalar superpartner masses are 
typically heavier than 1 TeV (compare table ^ but there is a reasonable chance for at least 
some squarks (but probably no sleptons) to be seen at the LHC. On the other hand, the 
gluino, while also preferably heavy, should have a much better chance of discovery at the 
LHC. Prospects for detection are also promising in direct detection searches for dark matter 
which are sensitive to spin-independent interactions. An improvement in sensitivity down 
to (Tp ^ ~ 10~^ pb will allow CDMS-II and other experiments to reach down to the bulk of 
the values favored by the posterior pdf (for both signs of and actually to fully explore 
all best-fit regions in the preferred case of ^ > 0. On the other hand, an improvement in 
sensitivity of at least 3 orders of magnitude will be required before favored ranges of cross 
sections for spin-dependent neutralino-proton interactions are tested by experiment. 

We stress that some of our findings do depend on our choice of the (usual) variables 
"iQ) f^i/2i ^0 and, especially, tan/3 as CMSSM parameters over which we scan, and on the 
choice of taking fiat priors on those variables. Of particular relevance is the upper bound 



rriQ < 4 TeV. Other choices are possible and should be examined. (See refs. |13, 53] and 
Note Added below.) The choice of CMSSM parameters and flat priors that we have made 
in the present analysis and the comparison with the mean-quality-of-flt statistics are meant 
to facilitate comparison with flxed-grid studies using similar assumptions. 

Note added: Very recently, after our analysis was completed, a new paper of Allanach 
et al., [^] has appeared. The authors argue that, instead of assuming a flat prior in tan/3, 
it is more natural to use a "REWSB prior" where ^ and the bi-linear soft mass parameter 
B are taken as inputs with flat distributions. Whether this choice (originally advocated 
by R. Ratazzi) is superior to any other is debatable but it is certainly justiflable to apply 
it, at least for the sake of examining the sensitivity of observables to the choice of priors. 
We note that with the REWSB prior the preference for large mo, well above 1 TeV, still 



remains. (A more detailed comparison is difficult because in ref. |53| a previous SM value of 
BR{B Xs^) has been used.) On the other hand, the authors in addition choose to give 
strong preference to cases where all the CMSSM mass parameters are of the same order. 
In our opinion this assumption does refiect a certain level of theoretical bias which at the 
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end strongly changes the conclusions obtained with the REWSB prior only. (In particular 
it disfavors the focus point region which, as we have shown, can be seen as being favored 
by the current results on BR{B — > Xg^).) We would prefer to see the lack of the hierarchy 
of the CMSSM mass parameters to be an outcome of applying experimental constraints, 
rather than of applying theoretical prejudice. 
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